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ABSTRACT 

We utilize the bulk viscosity of interacting strange quark matter to reevaluate the damping 
time scale. The presence of medium effect of bulk viscosity leads to a stronger damping of r- 
modes, which can be over an order of magnitude for realistic parameters. We find that the r-mode 
instability window is narrowed due to the medium effect, and hence when a pulsar reaches the 
instability window it will only slow down by gravitational wave emission to a period of 1.78msec 
instead of 2.5msec given by early estimate. As a theoretical upper rotation limit of pulsars, the 
period of 1.78msec is very close to the two most rapidly spinning pulsars known, with periods of 
about 1.6msec. 

PACS numbers: 97.60.Jd,12.38.Mh, 97.60. Gb 

Subject headings: dense matter-gravitation-stars:strange-stars:rotation 



1. INTRODUCTION 

At high density, normal hadron matter has been 
predicted a deconfined transition, where quark- 
gluon plasma is formed. This could have impor- 
tant consequences for compact stars, with central 
densities several times the nuclear saturation den- 
sity. Ever since strange quark matter (SQM), a 
conglomerate of up, down and strange quarks, has 
been suggested as a possible absolutely stable or 
metastable phase of nuclear matter(Witten 1984), 
it has been speculated that strange stars might ex- 
ist in the universe(Alcock et al. 1986; Haensel et 
al. 1986; Colpi and Miller 1992). 

The studies that quark interactions within low- 
est order perturbative QCD were considered in 
the MIT bag model also predicted that suffi- 
ciently heavy strangelets might be absolutely sta- 
ble(Fahari and Jaffe 1984). If the SQM hypothe- 
sis is correct then some (perhaps all) pulsars may 
be strange stars(Alcock et al. 1986; Haensel et al. 



1986; Colpi and Miller 1992). Thus it would have 
obvious implications for our understanding of pul- 
sars if strange stars could exist. It is an impor- 
tant subject that probing possible observational 
evidences of the existence of strange stars by as- 
trophysical investigations. 

Since self-bound of SQM due to strong interac- 
tion is very important, the mass-radius relation of 
the assumed strange stars has been uncovered to 
be different from that of neutron stars. However, 
for the canonical mass of 1.4M Q , gravity domi- 
nates the strong interaction which leads to strange 
stars and neutrons being similar in size(Alcock et 
al. 1986). Therefore, the attempt to distinguish 
strange stars from neutron stars for given observed 
masses and radii of pulsars seems to have encoun- 
tered with great difficulties. Strange stars could 
indicate the distinguishable signal from neutron 
stars by their cooling properties, but the difference 
in the cooling behavior of strange stars and neu- 
tron stars will disappears if the direct Urea process 



is operating in the core of an neutron star, which 
is done either by the proton fraction rising above 
10% (Lattimer ct al. 1994) or by hyperons being 
present in the core of a neutron star(Page et al. 
2000). This would make it difficult to identify a 
strange star from cooling data. 

Ever since Andersson realized that the r-modes 
are unstable at all rates of rotation in perfect 
fluid stars( Andersson 1998), a series of papers 
have investigated the many implications for grav- 
itational radiation detection and the evolution of 
pulsar (Friedman and Morsink 1998; Lindblom et 
al. 1998; Owen ct al. 1998; Kojima 1998; Madsen 
1998; Andersson et al. 1999a,b; Ho 2000; Reza- 
nia and Maartens 2000; Madsen 2000; Andersson 
et al. 2001). In recent years, several crucial is- 
sues regarding the astrophysical relevance of the 
r-mode instability have been investigated. Key re- 
sults concern the interaction between oscillations 
in core fluid and the crust (Bildsten & Ushomirsky 
2000; Andersson et al. 2000; Lindblom et al. 
2000), the role of the magnetic ficld(Spruit 1999; 
Rczzolla et al. 2000; Mendcll 2001) superfluid- 
ity(Lindblom & Mendell 2000; Andersson and 
Comer 2001) and the effect of exotic particles 
that are thought to exist in the deep neutron star 
core(Jones 2001; Lindblom et al. 2001). 

Meanwhile, Madsen (1998) has pointed out that 
the r-mode instability may provide the means to 
distinguish strange stars from neutron stars. The 
main reason for this is that the viscosity coef- 
ficients arc rather different in these two cases. 
While the shear viscosity of a strange star is com- 
parable to that of a neutron star, the bulk viscos- 
ity would be many orders of magnitude stronger 
than its neutron star counterpart. This has in- 
teresting effects on the r-modc instability. Based 
on the characteristic r-mode instability window of 
strange stars, which is related to the gravitational 
wave emission and the viscosity, a few mechanisms 
have been discussed to explain the clustering of 
spin-frequencies of LMXBs (Madsen 2000; Ander- 
sson et al. 2001). 

However, the bulk viscosity coefficient in Mad- 
sen's studies takes the one for the case of non- 
interacting quark gas in MIT bag(Madsen 1992). 
In fact, if the interaction among quarks is consid- 
ered, the bulk viscosity increases, which is calcu- 
lated by Zheng et al. (2002). As we will show in 
the following, the r-mode instability window for 



strange stars will be significantly modified, when 
the improved bulk viscosity is adopted. 

As demonstrated blow, we find the r-mode 
instability window is evidently narrowed due to 
the inclusion of medium effect (interactions among 
constituent particles). The medium effect makes 
the lowest critical spin frequency rise to 558Hz 
from about 300-400Hz(for the case of free quark 
gas). The limiting period of 1.78msec correspond- 
ing 558Hz is closer to the two most rapidly spin- 
ning pulsars known (with the periods of 1.56msec 
and 1.61msec) than the periods of 2.5msec and 
3msec given by Madsen. 

In addition, if it could be established that 
rapidly spinning pulsar are strange stars, the pul- 
sar data would put certain constraints on the 
model parameter, namely the current mass of 
strange quark mass m s . As we can see in Madsen's 
article (Madsen 2000), in the normal and especially 
in the 2-flavor color superconducting phase (2SC) 
case, the theoretical results is not consistent with 
the pulsar data when m s — lOOMeV. However, 
when we taking the medium effects into account, 
the constraints to m s from pulsar data appears to 
be deeply relaxed, and we will show that m s can 
be taken as lOOMeV in this article. 

This paper is organized as the following. In 
section 2, we recall the bulk viscosity of inter- 
acting quark matter and compute the damping 
time scale on r-modes. In section 3, we give 
the improved critical curve (spin frequency) in spin 
frcqucncy(i^)-temperature(T) plane. In section 4, 
we give a summary. 

2. BULK VISCOSITY AND DISSIPA- 
TION ON r-MODES 

Ever since Wang & Lu (1984) pointed out that 
strange quark matter is characterized by a huge 
bulk viscosity relative to nuclear matter, some in- 
vestigations have tried to calculate the relevant 
viscosity coefficient of SQM(Sawyer 1989; Madsen 
1992; Goyal et al. 1994). In a MIT bag model, it is 
thought that the exact solution had been obtained 
by Madsen(Madsen 1992). However, he ignored 
the coupling among quarks in the bag, which had 
been considered in the study of the equation of 
state of SQM(Schertler et al. 1997). Soon before, 
we found the coupling's effect on the bulk viscos- 
ity leads to an increase of the viscosity over an 
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order of magnitude although the medium modifi- 
cations of the equation of state of strange quark 
matter was proven to be small(Zhcng et al. 2002). 
To very good approximation, the relevant bulk vis- 
cosity coefficient still takes the form formulated by 
Madsen (1992): 
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but here a and (3 were given by (Zheng et al. 2002), 
which are extremely differing from Madsen's and 
strongly depending on the coupling constant of 
strong interactions among quarks. We expressed 
them as 
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where k Fi = {p 2 - mf ) 1/2 ,Ci = m - m*^- 
and to* was given in the references(Klimov 1982; 
Weldon 1993; Pisarski 1989; Blaizot and Ollitrault 
1993; Vija and Thoma 1995). 

The bulk viscosity is expected to be the dom- 
inated internal fluid mechanism in hot com- 
pact stars. The timescales for the bulk viscos- 
ity damping of r-modes need to be estimated. 
Since the coupling of the r-modes to bulk vis- 
cosity vanishes in the lowest-order expression, 
the derived timescales must be based on the 
fully self-consistent second-order calculation of 
this coupling. The earlier estimates can't gain 
this end(Lindblom et al. 1998; Andersson et al. 
1999a; Kokkotas and Stergioulas 1999). But it 
has been completed by Lindblom et al. (1999). 
We substitute our viscosity for the given Low-T 
limit(T < 10 9 K) into the formula(6.2)given by 
Lindblom et al. (1999), and find the only dif- 
ference comparing to the early calculation is the 
change of the viscosity coefficient. Thus, we can 
immediately obtain the timescale by a simple com- 
parison instead of the repeat of the complicated 
calculations 



t b = f B s(irGp/n 2 )T g - 2 , 
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with 



where Tg and toioo denote temperature in units of 
10 9 K and the current mass of strange quark in 
units of 100 MeV , and p is the mean density of 
the star. 

Evidently, f B is determined with the chemical 
potential p,d, which can be obtained by solving the 
equations related to chemical equilibrium, electric 
charge neutrality and conservation of baryon num- 
ber for given p, namely 
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f B = 2.83 x 10 3 a- 1 pm^ 100 . (5) 



Here, n denotes the baryon number density, and 
n-i = -^ik 2 Fi is the particle number density, is 
the energy per baryon and we approximately take 
it as the mass of a neutron in our calculation. 

The coupling g will extremely influences the so- 
lutions of the above equations. Following Schcrtler 
et al. (1997), g is taken as a free parameter rang- 
ing from to 5 and the equations(6) (7) (8) are 
solved numerically. Figure 1 shows f B as a func- 
tion of p for different g and m s — 200MeV. The 
time scales depend on g remarkably while they arc 
nearly independent of p. The time scale for g = 5 
arrives at tens of times shorter than the case of 
g = 0. For a strange star with mass M — 1.4M© 
and radius R = 10km, we find f B ranges from 
4.24xl0- 2 s( 3 = 0) to 1.53xl0~ 3 s( 3 = 5). This 
is of great interest because this will significantly 
increases the critical rotation angular velocity for 
the onset of r-mode instability. 

3. CRITICAL SPIN FREQUENCY FOR 
THE ONSET OF THE r-MODE IN- 
STABILITY 

The r-mode unstable (or stable) regime of the 
relativistic stars, neutron stars as well as strange 
stars, depends on the competition between the 
gravitational radiation and various dissipation 
mechanisms. To plot the instability window of 
r-mode or gain the critical rotation frequency for 
a given stellar model as a function of temperature, 
we need to acquire the characteristic timescales, 
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damping and growing timcscalcs of r-mode insta- 
bility. Due to the emission of gravitational waves, 
the r-mode grows on a timescale 

T G = T G (TTGp/n 2 ) 3 (10) 

where tq is — 3.26(1.57)s for n = 1(0) polytropic 
EOS, which was respectively studied by Lindblom 
et al. (1999) and Kokkotas and Stergioulas (1999). 
For the viscous damping timescales, we shall also 
consider the shear viscosity besides bulk viscosity 
discussed in the last section. In strange stars, the 
time scale for the shear viscous damping is given 

T5=f5(a s /0.1) 5/3 T 9 5/3 (11) 

Here, fs is 5.37(2.40) x 10 8 s corresponding to n = 
1(0), as is the strong coupling and we take as = 
0.1 in following calculation(Madsen 2000). 

We can now evaluate the critical spin frequency 
as a function of temperature from the equation 

— + — + — = (12) 

TG Ts TB 

Figure 2 shows the regions of r-mode (in)stability 
in spin frequency-temperature^ — T)planc for a 
strange star with mass M = 1.4M Q and radius 
R = 10km. The shading between the two curves 
displays the effect of the medium modification 
of quark masses on critical rotation frequencies. 
The medium effect narrows the r-mode instability 
window. The dotted curve corresponds to g = 0, 
reduced to the Madsen's result(Madsen 2000), and 
g = 5 is assumed in the upper contour(the solid 
curve). The upper contour has a lowest limiting 
frequency denoted by C as 558Hz(the correspond- 
ing period is 1.78msec), which is more close to the 
two most rapidly spinning pulsars known, with 
frequencies of 642Hz and 622Hz( the periods are 
1.56 and 1.61msec), than the period of 2.5msec. 
This implies a strange star would slow down by 
gravitational wave emission when it reaches the 
instability window and spin around in 1.78msec 
instead of the 2.5~3msec expected by Madsen. 
Figure 3 shows the results for 2SC stars. Simi- 
lar to figure 1, the medium effect also increases 
the critical rotation frequency of 2SC stars, but it 
should be stressed that the medium effect leaves 
the most rapidly stars away from the instability 
window. 



Figure 4 depicts the instability windows where 
the current mass of strange quark takes 100 MeV. 
If the rapidly rotation pulsars could be regarded 
as strange stars, the medium effect would relax the 
stringent constraint on the choice of QCD param- 
eters in contrast to the non-interacting medium 
case, for example, the 2SC stars for smaller cur- 
rent mass m s can safely exist due to the medium 
effect (see the upper dashed curve in figure 4). 

4. CONCLUSIONS AND DISCUSSIONS 

We apply the bulk viscous coefficient including 
medium effect to reevaluate and discuss the vis- 
cous damping time. The time greatly depends on 
the strong coupling g. We find the time scale tb, 
rather than non-interacting perfect fluid where fs 
is constant (The fact is that fs versus p slowly in- 
creases for non-vanishing m s as shown in figure 
1), weakly decreases with increasing mass density 
of stars for the given larger g. Therefore, the 
medium effect due to strong interactions among 
quarks adds the viscous dissipation of instability 
modes. For a star with M = 1.4M© and radius 
R = 10km, we numerically calculate the critical 
spin frequency as a function of temperature. The 
medium effect lead to a significant lift of the crit- 
ical frequencies. We find that the largest critical 
period is reduced to 1.78msec from 2.5msec, ap- 
pearing to be more close to the two most rapidly 
spinning pulsars known. 

We also study the 2SC stars and give our im- 
proved scenario. The current mass of strange 
quark in the model have a broader range taking 
the medium effects into account. 

Finally it should be mentioned that we here 
only consider weak coupling of strong interactions 
among quarks in MIT bag. The long range unper- 
turbation interactions are contained in the usual 
bag constant. Although the MIT bag constant is 
able to contribute to EOS of SQM, the long range 
effects on dynamical quantities, such as viscous co- 
efficient etc., seem to be ignored due to the hiding 
of the microscopic processes into the bag constant. 
If the processes can be taken into account, we may 
conjecture that the limiting spin of strange stars 
would shift further upward. This is our future 
work. 
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Fig. 1. — The time scales tb as functions of mass 
density for different couplings <?, which are 0, 1, 2, 
3, 4, 5 respectively from top to bottom 




Fig. 2. — Critical spin frequencies for strange 
stars as functions of temperature. The lower con- 
tour(dotted) curves stand for free quark gas. The 
shadow displays the medium effect. 
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Fig. 3. — As Fig 2, but assuming a 2SC state. 



400 



200 





' 1 ' 


1 — 


1 'f 1 


PI™ 


1H— r- 


I ' 1 










.'x/x 




























J • 


J 











s 


10" 


10 7 


10* 


10* 


10™ 


10 



T(K) 



Fig. 4. — The r-mode instability windows for the 
case of m s =100MeV. The solid curves denote the 
strange stars with normal phase and the dashed 
curves represent the case with 2SC phase. The 
lower solid and dashed curves correspond to the 
case of non-interacting quark matter. 
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